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ABSTRACT 

We analyze archived Chandra and XMM-Newton X-ray observations of 536 Sloan Digital Sky Survey 
(SDSS) Data Release 5 (DR5) quasars (QSOs) at 1.7 < z < 2.7 in order to characterize the relative 
UV and X-ray spectral properties of QSOs that do not have broad UV absorption lines (BALs). We 
constrain the fraction of X-ray weak, non-BAL QSOs and find that such objects are rare; for example, 
sources underluminous by a factor of 10 comprise <2% of optically-selected SDSS QSOs. X-ray 
luminosities vary with respect to UV emission by a factor of <2 over several years for most sources. UV 
continuum reddening and the presence of narrow-line absorbing systems are not strongly associated 
with X-ray weakness in our sample. X-ray brightness is significantly correlated with UV emission 
line properties, so that relatively X-ray weak, non-BAL QSOs generally have weaker, blueshifted 
C IV A1549 emission and broader C III] A1909 lines. The C IV emission line strength depends on both 
UV and X-ray luminosity, suggesting that the physical mechanism driving the global Baldwin effect 
is also associated with X-ray emission. 

Subject headings: galaxies: active — galaxies: nuclei — X-Rays: general — quasars: absorption lines 
— quasars: emission lines 



1. INTRODUCTION 

It has been known for some time that the UV and 
X-ra y luminosities of quasars (QSOs) are correlated 
(e.g. JAvni &: TananbaunJ|19 82'. and references therein), 
and recent studies have carefully quantified this relation 
across «5 orders of magnitude in UV luminosity (e.g ., 
IStrateva et al]|2005t ISteffen et al1l2006l : lJust et alboOTl ) . 
Such studies inform ongoing efforts to understand the 
structure and physics of QSO nuclear regions, provid- 
ing quantitative constraints on models of physical asso- 
ciations between UV and X-ray emission. Because UV 
photons are generally believed to be radiated from the 
QSO accretion di sk while X-rays are prod uced in the 
disk corona (e.g., iRevnolds fc NowaS l2003l and refer- 
ences therein), the UV/X-ray luminosity relation is an 
indication of the balance between accretion disks and 
their coronae. For example, a large fraction of intrin- 
sically X-ray weak sources would suggest that coronae 
may frequently be absent or disrupted in QSOs. In this 
study, we combine results from recent optical/UV and 
X-ray observations of hundreds of QSOs in order to con- 
strain the fraction of sources that are anomalously X-ray 
weak, and also to test for additional physical effects that 
contribute to the scatter observed in the relation between 
UV and X-ray luminosities. 

Radio-lou d QSOs are well known to be relatively X-ra y 
bright (e.g.. lWorrall et al.lll987l : iBrinkmann et al.ll2000D . 
while QSOs with broad UV absorp t ion lines (BALs) ar e 
X-ray faint fe.g.. lGreen et al.lll995t fBrandt et al.l l200(/). 
We wish to understand the emission processes of QSOs 
without the added complexity of the strong X-ray ab- 
sorption associated with UV BAL outflows or the en- 
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hanced X-ray emission that may be associated with a ra- 
dio jet. In order to accomplish this, BAL and radio-loud 
QSOs must be carefully removed from samples. This pro- 
cess can be complicated in cases where available spectra 
do not extend to wavelengths where strong UV BAL ab- 
sorption may occur. 

In this study, we compare the UV and X-ray properties 
of optically-sele cted, radio-quiet Sloan Digital Sky Sur- 
vey (SDSS; e.g.. lYork et"alll2000( l QSOs that are known 
not to host BAL outflows along the line of sight. Careful 
screening enables us to study the ways in which these 
"ordinary" QSOs deviate from the general relation be- 
tween UV and X-ray luminosities. We quantify the scat- 
ter about the best-flt relation between UV and X-ray 
luminosities, and in particular constrain the fraction of 
optically-selected QSOs that are intrinsically weak X-ray 
emitters. Secondly, we search for correlations between 
UV emission/absorption properties and relative X-ray 
brightness, where the term "relative" indicates the ob- 
served X-ray brightness compared to that expected for 
an average QSO with the same UV luminosity. Beyond 
the UV/X-ray luminosity relation, our flndings place ad- 
ditional constraints on physical models that relate QSO 
UV and X-ray emission processes. 

The significant amount of scatter in the X-ray bright- 
ness of individual sources compared to that of "average" 
QSO s with the same UV luminosity (e.g., Steffen et aT] 
|2006[) indicates that some unmodeled physical factors in- 
fluence the relation between UV and X-ray luminosities. 
Besides the cases of radio-loud and BAL QSOs men- 
tioned above, possible additional causes of scatter in- 
clude additional X-ray absorption that is not associated 
with UV BALs, intrinsically weak (or strong) X-ray emis- 
sion, and time variability. Studies of X-ray brightness in 
QSOs use different methods to quantify UV absorption 
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and X-ray weakness (complicating comparisons between 
different samples), and intensive UV and X-ray spectro- 
scopic campaigns have been performed for only a few 
QSOs. As a result, we do not have strong constraints 
on how frequently X-ray absorption, intrinsically faint 
X-ray emission, and time variability cause X-ray weak- 
ness in QSO samples. In the remainder of this section, 
we briefly review published cases involving these physi- 
cal scenarios. This discussion inclu des all the non-BAL 
QSOs listed in iBrandt et al .' (2000) and five of the six 
sources from I Yuan et al.l ( 1998,1 with aox < —2.^ We 
also discuss cases of unusual X-ray weakness in recent 
studies of individual sources. 

1.1. X-Ray Absorption 

Even for QSOs that do not show UV BAL absorp- 
tion, strong X-ray absorption may be the simplest expla- 
nation for X-ray weakness. Mrk 304 (PC 2214 + 139) 
is an X-ray weak source that is often classified as a 
Seyfert 1 gala xy. It does not show UV BAL absorp- 
tion, although iLaor fc Bra^ (|2002[ ) found some evi- 
dence for moderate C IV and N V line absorption. In the 
X-rays, it is strongly absorbed by a multizone, ionized 
absorber with a total column density Nh ^ 10^^ cm~^ 
(jPiconceUi et al.ll2004l : [Brinkmann et al.ll2004D . Another 
X-ray weak source, PC 1126 — 041, shows evidence for 
an ionized X-ray absorber. However, low- velocity, broad 
C IV absorption has been found in some observations, 
indicating that PG 1 126 — 041 should b e classified as a 
(variable) BAL QSO (jWang et al.lfl999l ). 

1.2. Intrinsic X-Ray Weakness 

X-ray weakness cannot always be attributed to absorp- 
tion. The nearby (z = 0.192), narrow-line type 1 QSO 
PHL 1811 shows no evidence of BALs in its UV spec- 
trum. With an upper limit for a neutral X-ray absorb- 
ing column of Nh < 8.7 x 10^" cm~^, the sour ce is not 
strongly X-ray absorbed (|Leighlv et al.l [200 7bf) . It has 
been consistently X-ray weak since it was first observed 
almost 20 years ago with ROSAT, although the X-ray 
flux has varied by a factor w5. This variation (small 
compared to the degr ee of X-ray weakness) suggested to 
iLeighlv et all (|2007b[ ) that the X -ray emission was not 
scattered into view, as the medium responsible for the 
scattering would need to be implausibly small. PHL 1811 
therefore appears to be an intrinsically weak X-ray emit- 
ter. T he optical/UV spectr um of PHL 1811 is also un- 
usual (|Leighlv et al.ll2007a[ ). It is very blue and shows 
no forbidden or semi-forbidden line emission. The C IV 
A1549 line emission is weak by a factor of ~ 5 comp ared 
to the composite spectrum of iFrancis et al.l ()1991|) . In 
this study, we will test whether intrinsically X-ray weak 
objects like PHL 1811 (which would have been flagged 
for SDSS spectroscopy based on its optical colors alone) 
are common, and whether these unusual emission line 
characteristics are associated with X-ray weakness. 

Another source, P G 1011—040, is relatively X-ray weak 
by a factor of ~10 (jBrandt et al.l 120001 : iGallagher erall 
|2001|) . There is some evidence for C IV absorption in 

^ The si xth source, QSO 03 16 — 346, was identified from optical 
spectra bv lMonk et al.l II1988I) . UV spectra are required to search 
for BAL absorption, but no UV spectra of this source are present 
in the MAST or HEASARC archives. 



the UV spectrum, but no strong indicati on of a BAL 
(|Brandt et al.ll2000l : iLaor fc Brandtl l2002l ) . It appears 
to be relatively unabsorbed in X-rays, with an upper 
limit on the absorbing column density of Nh < 5 x 
10^^ cm^^. In contrast to PHL 1811, multiple obser- 
vations of PG 1011 — 040 have not found much X-ray 
luminosity variation (|Gallagher et al.ll2001h . 

1.3. UV and X-Ray Variability 

Source variability may be an important factor con- 
tributing to the scatter in the UV/X-ray luminosity rela- 
tion. Additional scatter is introduced by the fact that op- 
tical/UV and X-ray observations are usually taken some 
time apart, and this time may reach up to years. 

The narrow-line Seyfert 1 galaxy Mrk 335 has recently 
been reported to have decreased in X-ray brightness by 
a factor of «30 (jGrupe et al.ll2007ah . The authors sug- 
gested the variation may have been caused by the onset 
of heavy X-ray absorption, which may also lead to the 
appearance of UV BALs. PG 0844 -f 349 has been ob- 
served to vary in the X-rays by up to 60% on short (20 ks) 
time scales, and by a factor of ^10 on longer (multi-year) 
time scales (e.g.. IGallagher et al.ll200l] : iBrinkmann et al.l 
120031 ). The narrow-line Seyfert 1 galaxy WPVS 007 has 
dimmed in the X-rays by a factor of 100 or more, while 
a C IV BAL a ppears to be formin g in the UV spectrum 
of this source (jGrupe et al1l2007bl ^. 

Finally, we note that our understanding of the causes 
of X-ray weakness is limited by observing time scales. 
Over longer time scales. X-ray weak QSOs may show dif- 
ferent properties (such as long-term emission or absorp- 
tion variation), and may even cease to be X-ray weak. 
We have also noted in §?? cases where UV BALs may 
be transient, and recent study has demonstrated that 
BALs evolv e on multi-year (rest- frame) time scales (e.g., 
iGibson et~al . 2008.) . For these reasons, we are only able 
to place limits on the extent of UV/X-ray variability 
( §3.5p . Future near-simultaneous, multi- wavelength ob- 
servations of QSOs will improve on these constraints. 

Throughout this work we use a cosmology in which 
Ho = 70 km s-i Mpc^^, = 0.3, and f^A = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

The SDSS Data Release 5 (DR5; 

lAdelman-McCarthv et all l2007l ) has recently cata- 
loged spectra of over 77,000 QSOs, most of which were 
selected for spectrosco py based on optical considerations 
([Schneider et al.ll2007| ). Many of these sources have also 
been observed in the X-rays with modern telescopes 
such as Chandra and XMM-Newton, enabling detailed 
studies of relations between UV and X-ray properties 
for a large sample of optically-selected QSOs. 

In this study, we consider the 536 SDSS QSOs at red- 
shifts 1.7 < z < 2.7 that have been observed with Chan- 
dra or XMM-Newton. The lower limit on redshift ensures 
that the SDSS spectrum reaches to «1400 A, allowing 
good coverage of any potential high-velocity C IV ab- 
sorption, up to « —30,000 km s~^. The upper redshift 
limit ensures that the 2500 A continuum is in view so 
that we can determine the UV continuum flux accurately. 
We o btain QSO redshifts fr om the SDSS DR5 QSO cat- 
alog ( ISchneider et al.l[2007[ ). These redshifts were gen- 
erated by the SDSS spectroscopic pipeline based on the 
overall continuum and emission properties of each spec- 
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trum; a small fraction were revised during vi- 

sual inspection for the QSO catalog. The mean error in 
redshift estimated from the pipeline for our sources is 
RiO.002. An rms difference of 0.006 in redshift determi- 
nation has been ob served in sources o bserved m ultiple 
times in the SDSS (jWilhite et al.ll2005t fSchneider et all 
l2007f) . We therefore expect that errors in redshift are 
typically 0.5% or less for our sources. We note recent 
studies have found no significant evidenc e for redshift 
evolution in an y fe.g.. ISteff en et al."2006l. X-ray spec- 
tral shape (e.g., IShemmer et al. 2005), or infrared con- 
tinuum properties (e.g. JJiang et al. 2006), so we do not 
expect that our results are strongly influenced by our 
redshift restriction. 

In order to illustrate the parameter space covered by 
our sample, we plot in Figure [1] the redshifts and ab- 
solute i magnitudes for our sources together with the 
full DR5 QSO catalog. For comparison, we also indicate 
SDSS DR5 QSOs at z < 0.5 which were included in the 
Bright Quasar Survey (BQS; ISchmidt fc Greenlll983f l in 
the plot. 

2.1. X-Ray Observation Catalog 

We obtained the list of all Chandra observations that 
were publicly available as of 2007 July 26 and searched 
for cases where a DR5 QSO fell within 15' of the tele- 
scope pointing for an ACIS observation with no grat- 
ings. For this subset of Chandra observations, we tested 
whether each QSO was on an ACIS chip using the 
CIAO tool dmcoords. We retrieved observations with on- 
chip sources from the Chandra online database and ran 
acis_process_events on each data set in order to apply 
the latest calibration data, including the time-dependent 
gain correction and charge transfer inefficiency (CTI) 
corrections.'^ All processing was performed using CIAO 
version 3.4.'^ 

We extracted spectra for the source and the back- 
ground of each on-chip DR5 QSO using the CIAO tool 
psextract. For the source extraction region, we used 
a circle with radius equal to the 90% encircled energy 
fraction plus an additional 5 pixels. For the background 
region, we used an annulus around the source with in- 
ner (outer) radii equal to the 90% encircled energy frac- 
tion plus 15 (50) pixels. We chose these regions because 
we have found that they accurately characterize the sur- 
rounding X-ray background with a minimum of contami- 
nation by nearby sources in the background regions. Be- 
cause a large fraction of our sources are serendipitously 
observed at large off-axis angles and may be faint in any 
case, we are concerned to avoid this contamination. 

Because the CIAO tool psextract does not generate 
spectra for sources with zero counts, we flagged these 
cases for proper handling by our fitting algorithm de- 
scribed in i|2.3l We also excluded observations where the 
source lies within 32 pix els of a chip edge. 

We used wavdetect (jFreeman et al.l [2002( 1 with the 
sigthresh parameter set to 10~^ to detect any sources 
that may contaminate the annular regions we used to 
extract background spectra. In six cases, we found that 
the total number of counts contributed to the background 



^ See t he CIAO thre ads|http://a sc. harvard.edu /ciao/threads/acistim ei 
|http: //asc. ha rvard . edu/ ciab / threads / acisapplycti/] 
http:/ /asc. harvard.edu/ciao3. 4/ 



spectrum by detected sources was > 10% of the observed 
spectrum. In these cases, we manually adjusted the back- 
ground extraction regions to be free of source contami- 
nation and re-processed the data. 

We also obtained a list of all XMM-Newton observa- 
tions that were public as of 2007 August 22 for which a 
DR5 QSO lay within 15' of the target. We reduced EPIC 
camera data using the emchain and epchain tasks in the 
SAS 7.1.0 package.** We removed intervals of background 
flaring from each observation by flagging all time inter- 
vals for which the count rate exceeded the baseline count 
rate by a factor of three or more, where the baseline count 
rate was defined as the rate for which the distribution of 
count rates peaked. 

We extracted spectra for the MOS and pn cameras 
using a circular aperture 15" in radius for the source 
and an annulus with inner (outer) radii of 50" (75"). 
For each arcminute that the source was off-axis from the 
nominal pointing, we increased the extraction radii by 
1" (2.3") for the MOS {pn) cameras to account for the 
increase in the constant encircled energy fraction radius. 
As for the Chandra processing, the extraction regions 
were chosen to characterize effectively the surrounding 
background while minimizing contamination from nearby 
sources. We used the SAS utility eboxdetect to detect 
sources that may contaminate our background regions. 
In cases where bright sources were detected in our back- 
ground extraction regions, we omitted the the brightest 
source from the background extraction in order to im- 
prove background estimation. 

A single XMM-Newton observation may include one or 
more exposures from the MOSl, M0S2, and pn cameras. 
From the available set of exposures, we preferentially se- 
lected the MOS observation with the longest exposure. 
Although the pn cameras provide a larger effective area, 
we found that the higher angular resolution of the MOS 
cameras yield more reliable background estimation and 
a higher fraction of source detections. For this reason, 
we prefer the MOS cameras for our study. 

We discard sources for which the effective area is 
<100 cm'^ in all bins in order to filter out sources that are 
near a chip edge. Even so, we find that the fraction of 
undetected sources is much greater in the XMM-Newton 
reductions than for Chandra, due to the significantly bet- 
ter angular resolution of the Chandra telescope. In cases 
where undetected sources are important for our statistics, 
we exclude the XMM-Newton observations from consid- 
eration. 

In eight cases, the number of background counts at- 
tributed to eboxdetect sources in the background ex- 
traction regions was >30% of the total counts extracted 
from the source region. Visual inspection indicates that 
most of these are cases where the source was weak and 
the candidate background sources are borderline detec- 
tions. We manually selected a suitable background and 
re-processed the two cases where strong sources clearly 
impact the background estimation. 

Approximately 9% of our Chandra QSOs and 1% of our 
XMM-Newton QSOs were targets of observations. (The 
fraction is higher for Chandra because targets observed 
with both telescopes will only have their Chandra ob- 
servations included in our sample.) While there may be 

jam/ and r- / 

■* http://xmm.vilspa.csa.es/sas/ 
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small biases introduced by the tendency of observers to 
select unusual objects for observation, our study is dom- 
inated by sources observed serendipitously. We remove 
targeted sources from our sample when appropriate. 

In all, we have X-ray observations of 536 SDSS QSOs 
at 1.7 < z < 2.7. Of these, 315 were obtained with Chan- 
dra, while 221 were obtained with the XMM-Newton 
MOS cameras. We considered any source that was de- 
tected (using Poisson statistics) at 99% confidence above 
the background in the soft (0.5-2 keV), hard (2-8 keV), 
or full (0.5-8 keV) bands to be detected in X-rays. 90% 
of the Chandra sources and 67% of the XMM-Newton 
sources in our sample were detected by this criterion. We 
provide a catalog of our X-ray observations in Table [TJ 

2.2. Fitting SDSS Spectra 

Before fitting the SDSS spectra, we multiply them 
by a constant to match (approximately) the r, and 
i PSF magnitudes synthesized from the spectra to the 
measured photometric magnitudes. We also correct the 
sp ectra for Galac tic extinction using the reddening curve 
oflCardelli et al.l ([1989) with the near-UV extension of 
lO'DonneU (11991 . Wc obtain E(B - V) from the NASA 
Extraga lactic Datab ase J|NED)^, which uses the dust 
maps of lSchlegel eTalJ (l998( ). The QSOs in our study 



are at such high redshifts and luminosities that their 
spectra are not significantly contaminated by their host 
galaxies. 

We fit SDSS spectra using the algorithm of 
iGibson et al] (|2008f l. which we summarize here. Our con- 
tinuum model is a power law redden ed using t he Small 
Magellanic Cloud reddening curve of m (fl99l . We do 
not attach physical significance to the intrinsic E{B — V) 
values we obtain from these fits, as the reddening is de- 
generate with the shape of the underlying continuum 
emission. The UV luminosities we report are therefore 
corrected for Galactic, but not intrinsic, reddening. This 
follows the general practice for studies of UV/X-ray lu- 
minosity relations. We fit regions that are generally free 
from strong absorption or emission features (1250-1350, 
1600-1800, 1950-2050, 2150-2250, and 2950-3700 A), 
and then re-fit the spectrum iteratively. At each itera- 
tion, we ignore wavelength bins that deviate by >3a from 
the continuum in order to account for unmodeled absorp- 
tion and emission features. We then fit Voigt profiles to 
the strongest emission lines expected in the spectrum: 
Si IVA1400, C IVA1549, C III]A1909, and Mg IIA2799. 
Our line wavelengths are taken from the SDSS vacuum 
wavelength list.® We fit emission lines iteratively as well, 
ignoring at each step bins that are absorbed by more 
than 2.5(7 from the continuum -I- emission fit. 

We searched for narrow absorption line (NAL) features 
attributable to the Mg II A2799 doublet in each spec- 
trum. We searched each spectrum from 1550-2799 A and 
fiagged a candidate Mg II NAL feature in cases where 
spectral bins were found at the expected doublet separa- 
tion which lay at least Aa and 3a below the continuum 
for the red and blue sides of the doublet, respectively. In 
order to screen out broader absorption features, we also 
required that the spectrum reached a level within 3.5(7 of 
the continuum in the bin centered between the red and 

^ http:/ /nedwww. ipac.caltech.edu/ 

® http:/ /www. sdss.org/dr6/algorithms/linestable. html 



blue components. We also required that the equivalent 
width (EW) of at least one of the red and blue compo- 
nents be < —0.5 A, in order to ensure that the features 
were reliably detected. Our results are not strongly de- 
pendent on the particular value of the EW threshold. 

We searched for C IV BALs by calculating BIq, which 
is similar to th e trad itional "balnicity index," BI, of 
IWevmann et al.l (|199lD . except that it is integrated from 
velocity offsets -25,000 km s~^ to km from the 
QSO rest frame. BI is roughly a measure of the EW 
of features that are absorbed across a span of at least 
2000 km s~^. Sources with BIq > are considered to 
have BALs. We tested for Mg II A 2799 BALs in a similar 
fashion (when spectra extended to 2800 A) , although the 
continuum in that region is difficult to determine due to 
Fe II emission and structure in the Mg II emission line. 
We found no cases where Mg II BALs were present with- 
out C IV BALs, and therefore use BIq measured for the 
C IV absorption region to distinguish between BAL and 
non-BAL QSOs. 

The UV properties of our sources are listed in Table [1] 
and Table El 

2.3. Fitting X-Ray Spectra 

We fit X-ray spectra for each source in order to find 
and constrain the monochromatic luminosity at 2 keV, 
L2 kcV- To accomplish this, we fit a broken power law 
model to each spectrum, with the power law break fixed 
at 2 keV in the rest frame. This model accomodates 
possible soft X-ray absorption or excesses compared to 
the hard X-ray spectrum. In cases where the source is 
not much brighter than the background, we fix F = 2 for 
the power law fit to the entire X-ray band. We do this in 
cases that have fewer than 10 counts in the total (source 
plus background) spectrum, and also in cases that have 
no source spectral bins exceeding the background by at 
least a factor of two when the spectrum is adaptively 
binned to contain 10 counts per bin. Wc fixed F = 2 for 
109 sources in our full sample A, including 12 of the 139 
radio-quiet, non-BAL QSOs in sample B (described in 
§??). We do not use the values of F in our study, as they 
can often only be weakly constrained even when they are 
not fixed to F = 2. 

In all cases , we fit the unbinned data using the Cash 
statistic, C (jCashl Il979f) . For XMM-Newton observa- 
tions, we fit only the spectrum from the camera selected 
as described in We estimate the 68% confidence 

error in L2 keV by adjusting the model normalization pa- 
rameter upward and downward, re-fitting the power law 
photon indices (unless they are fixed to 2 for background- 
dominated spectra) at each step. We then estimate the 
upper and lower limits of L2 keV from the model fits at 
the points where AC = 1. 

Values of L2 kcV for our sources are listed in Table [TJ 

2.4. Estimating Radio Loudness 

All but 25 o f our 536 sources we re covered in the FIRST 
radio survey (jBecker et al.lll995[ ). For these sources, we 
obtain the 1.4 GHz core flux densities from the DR5 
QSO catalog. For sources that were not covered by the 
FIRST survey, we used 1.4 GHz fiux densi ties obtained 
from the NVSS survey (jCondon et al.|[T998i) . We esti- 
mate the monochromatic luminosity at 5 GHz assuming 
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the radio flux follows a power law with a spectral index 
a — —0.8. We then calculate the rad i o-loudness param- 
eter (e.g., ISramek &: WeedmanI 119801 : iKellermann et al.1 

[Mi), 

We classify sources with log(i?*) > 1 as "radio- loud" 
and all other sources as "radio-quiet." In 457 of 536 
cases, the radio surveys (FIRST or NVSS) were sensitive 
enough to measure radio luminosities down to log(i?*) = 
1. For 508 cases, the radio data are sensitive down to 
log(i?*) = 1.1, and in the worst case, the data are only 
sensitive to log(i?*) = 1.5. Of the 457 sources with radio 
observations sensitive down to log(i?*) = 1, 33 are radio- 
loud (with log(i?*) > 1), and 21 of these have log(i?*) > 
1.5. We therefore estimate that only a small number 
(«2) of radio-loud sources have contaminated our "radio- 
quiet" sample, and even these have log(i?*) < 1.5 and 
therefore are not highly radio-loud. 

3. ANALYSIS OF PHYSICAL PROPERTIES 

We have obtained Chandra or XMM-Newton X-ray 
spectra for 536 SDSS DR5 QSOs at redshifts 1.7 < z < 
2.7. In cases where the source was observed with both 
Chandra and XMM-Newton^ we use the Chandra obser- 
vation with the longest effective exposure. We call this 
set of sources sample A. We use sample A in cases where 
we are concerned with the general properties of a large 
number of sources. 

Many sources in sample A were not detected in X-rays. 
In order to increase the fraction of detected sources, 
we construct a subsample, which we call sample -B, of 
sources observed with Chandra that are <10' off-axis 
and have at least 2.5 ks of exposure. We do not in- 
clude XMM-Newton sources in sample B because the 
lower sensitivity leads to a larger fraction of X-ray non- 
detections. We also limit sample B to contain only 
sources from the SDSS DR5 QSO catalog that were se- 
lected for SDSS spectroscopy based on optical criteria. 
These sources include those flagged as "low-x," "high-2," 
"star," or "galaxy" accor ding to the ta rget selection de- 
scribe d in [Richards et al.l (j2002B l and IStoughton et al.l 
(l200l . Those sources that were included in the cata- 
log purely due to their radio or X-ray properties are ex- 
cluded from sample B. Finally, we exclude sources that 
lie within 1' of the Chandra target in order to remove bias 
due to observers' tendency to point at unusual sources. 
(Dropping this last criterion does not materially affect 
our results below.) Sample B contains 163 sources, 139 
of which are radio-quiet, non-BAL QSOs. All sources in 
sample B were detected by Chandra in X-rays. Table [3] 
summarizes the number and types of objects in samples 
A and B. 

Sample B is defined without direct reference to source 
properties and is therefore representative of optically- 
selected SDSS QSOs. We use sample A for tests of cor- 
relations involving X-ray properties, but use sample B in 
cases where non-detections adversely affect our statistics. 

3.1. Distribution of Aaox 
From our optical/UV and X-ray fits, we calculate the 
monochromatic luminosities ^osooA -^2 kcV a-t 2500 A 



and 2 keV, respectively. We adopt the traditional defini- 
tion: 

aox^ 0.3838 logfl^y (2) 

aox is a logarithmic measure of the X-ray brightness of 
a source relative to its UV brightness. For an ensemble 
of QSOs, aox is known to decrease as a function of QSO 
U V luminosit y . Thi s trend was most recently quantified 
bv lJust et all ()2007t ). who found 

^ox{L,,ooa) = (-0-140 ± 0.007) \og{L^,,,^) + (2.705 ± 0.21^) 

We use the EM (estimate and m aximize) regres- 
sion algorithm ([Dempster et al.l 119771 ) in the Astron- 
omy Survival Analysis (ASURV) software pa ckage (e.g., 
Ilsoije fc FeigelsonI llQQOl: iLavallev et aTi I1992D to fit the 
trend of aox with UV luminosity for radio-quiet, non- 
BAL QSOs in sample B. Since all such sources have 
X-ray detections in sample B, the E M algorithm gener- 
ates a traditional "least-squares" fit (jlsobe et al.lll98i) 
in this case. We find a relation somewhat steeper than 
Equation [3] 

«ox(i2500i) = (-0.217 ± 0.036) logiL^,,,^) + (5.075 ± l.llp|) 

The values of aox iov our sample, together with the fits 
from Equations [3| and [4l arc shown in Figure [2| 
We then define 

Aaox = aox - aox (^25001)' (5) 

where aox is the value measured from observations. 
Aaox quantifies the relative X-ray brightness of a source 
compared to that of "ordinary" QSOs. For example, 
Aaox = —0.5 corresponds to X-ray weakness by a factor 
of «20 compared to QSOs with the same UV luminosity. 

The median of the magnitude of difference in aox es- 
timated from UV luminosity using Equations [3| and [4] 
is 0.02 for our detected sources, which is less than the 
median measurement error in aox and is small com- 
pa red to the range of Aaox in our sample. The study 
of I Just et al.l p007) used a different sample of objects 
and different methods for estimating luminosities, so the 
small difference in results is not surprising. We use Equa- 
tion [4] throughout to determine Aaox in order to cor- 
rect for any small biases that would cause our results 
to differ. We find no evidence of a correlation between 
Aaox and z for radio-quiet, non-BAL QSOs that might 
arise from systematic biases in our data. Because our 
data cover a relatively small redshift range, this result 
does not strongly constrain redshift evolution for QSOs 
in general. 

Figure [3] shows the distribution of Aaox for our 
sources. We divide the sources into four categories, 
based on the presence (or absence) of a UV BAL and 
whether they are radio-quiet (log(i?*) < 1) or radio- 
loud (log(i?*) > 1). Sample A contains 433 radio-quiet, 
non-BAL QSOs, 34 radio- loud non-BAL QSOs, 64 radio- 
quiet BAL QSOs, and 5 radio-loud BAL QSOs. The 
second panel of Figure [3| demonstrates that radio-loud 
sources are relatively X-ray bright, while the third and 
fourth panels show that BAL QSOs are X-ray weak com- 
pared to non-BAL QSOs. The distribution of Aaox f^or 
radio-quiet, non-BAL QSOs in sample B is shown in Fig- 
ure m An Anderson-Darling test rejects the hypothesis 
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that Aaojf is normally distributed for radio-quiet, non- 
BAL sources in B at 97.5% confidence.'' 

Because the Anderson-Darling test has only indicated 
non-normality at 97.5% confidence, we cannot rule out 
the possibility that the und erlying di stribution is normal. 
In fact, we note that .Strateva et al.l ((2005.) found that a 
Gaussian fit provided a reasonable approximation to the 
distribution in their larger sample that covered a wider 
range of UV luminosities and was also slightly contami- 
nated by BAL QS Qs. However, we were not able to test 
quantitatively the iStrateva et al.l (|2005D distribution for 
normality due to the X-ray non-detections in their sam- 
ple. 

In order to investigate any possible non-normality of 
the distribution of Aaox for sample B, we have fit the 
distribution with a Gaussian using the IDL GAUSSFIT 
routine. The resulting fit, with mean fj, « 0.01 and 
a « 0.10, is shown in Figure ID Compared to this fit, 
the measured distribution perhaps has a broader base 
and a narrow core at Aapx ~ 0. The breadth at 
the base of the distribution may be due to unmodeled 
X-ray absorption in non-BAL sources in the X-ray weak 
wing. We also speculate that some sources in the X-ray 
bright wing which we did not consider formally radio- 
loud may nonetheless have enhanced X-ray emission as- 
sociated with a modest radio jet. Detailed study of indi- 
vidual sources is needed to determine the significance of 
these effects. 



3.2. How Common are Optically-Selected, X-ray Weak 
QSOs (XWQs)? 

We would like to understand the distribution of rela- 
tive X-ray brightness/ weakness in order to determine the 
frequency of intrinsically X-ray weak objects. If strong, 
intrinsic X-ray emission is not a universal property of 
QSOs, this would have significant implications for our 
understanding of AGN physics. The completeness of 
X-ray AGN surveys could also be negatively impacted 
by a significant fraction of intrinsically weak sources. 

Because we are able to identify BAL QSOs unam- 
biguously, our sample is well-suited to provide tight 
constraints on the distribution of relative X-ray bright- 
ness/weakness. All of the 139 optically-selected, radio- 
quiet, non-BAL sources in sample B are detected in 
X-rays according to the criteria in iJ2.1l 

The minimum Aaox in sample B is —0.37, corre- 
sponding to relative X-ray weakness by a factor of «9. If 
the probability of finding a QSO with Aaox < —0.37 is 
«1.6%, then there is a >90% chance we would observe at 
least one such QSOs in 139 attempts, using a binomial 
distribution. Since we have observed no such sources, we 
conclude that <1.6% of optically-selected, radio-quiet, 
non-BAL SDSS QSOs have Aaox < -0.37. 

We have also increased our sample of optically-selected 
QSOs to include the 87 Palomar -Green (PG ) QSOs from 
the Bright Quasar Survey fBOS: [Sch^^Iidrfc Greenlll983^ 
at z < 0.5. These sources are well-studied in the optical, 
UV, and X-ray bands. Figure [T] shows that the BQS 

''' The Anderson-Darling test, related to the Kolmogorov- 
Smirnov test, calculates a statistical probability that the observed 
data points could co me from a norm al distribution. We used case 
three of the test ('e.g.. lStephenslll974l ). in which the mean and width 
of the (putative) normal distribution are not known a priori. 



sources extend our coverage to lower luminosities as well 
as lower redshifts. 

We have obtained UV and X-ray monoc hromatic lu- 
minos ities for BQS sources from Table 2 of Steffe n et al.l 
(2006"). We exclude 8 sources known to host BALs. Five 
sources (0043 + 039, 1001 -f 054, 1004 + 130, 1700 + 518, 
an d 2112 + 059 ) have been identified as BAL QSOs 
bv lBrandt et al.l ()200G[ ) and the authors th ey reference. 
1351 -\- 640 also hosts a C IV BAL (e.g., IZheng et al ^ 



2001 1 ) . 1126-041 hosts a variable C IV BAL (jWang et~ 
19991 see esp. their Figure 1 ). We exclude 1411 -I- 442 



Ce.g.- lMalkan et a .1119871 : iBrandt et ani2000( l and 1535- 
547 (jSmith et al.lll997f ). both of which have broad ab- 
sorption at low velocities; the latter is also highly po- 
larized, similar to many BAL QSOs. For 1259 -I- 593, 
we use the X-ray monochromat ic luminosi ty obtained 
with XMM-Newton measured bv lFang etldl l|2008lV this 
source was previously undetected in less-sensitive X-ray 
observations. It is possible that our remaining sample of 
BQS QSOs is contaminated with a small number of BAL 
QSOs. Brandt et al. (2000) tested for C IV BAL absorp- 
tion in the spectra 55 of the 87 BQS QSOs at z < 0.5. 
The X-ray weak QSOs (with aox < -1.8) all had C IV 
coverage, an d in any case BAL QSO s make up a minority 
(«15%; e.g.. iHewett fc Foltz|[2Q0l of optically-selected 
samples, so any any remaining mis-classified BAL QSOs 
would have only a small influence on our statistics. 

There are three BQS sources with Aaox < —0.4, and 
the presence of these sources in our combined sample 
weakens our constraints on the fraction of X-ray weak 
QSOs. The sources are: PG 0844 -f 349 (which is highly 
X-ray variable), PG 1011 — 040 (which is relatively unab- 
sorbed in X-rays), and PG 2214-1- 139 (which is strongly 
X-ray absorbed). All three of these sources were dis- 
cussed in §?? as examples of unusual X-ray weakness. 

With these BQS sources included in our combined 
sample, we now have 3 sources (out of 217) with 
Aaox < —0.4, corresponding to a constraint that <3.1% 
of optically-selected QSOs have Aaox < —0.4. It is 
not clear why the BQS sources have a greater fraction 
of X-ray weaker QSOs than wc find in our pure SDSS 
sample. The cause may be random fluctuations in the 
data, or it may be due to physical differences in the 
sources, as the nearby BQS QSOs have lower UV lu- 
minosities than the SDSS QSOs in our sample. In any 
case, the differences in the upper limits we are able to 
draw are not great. The large number of sources in 
the combined SDSS and BQS sample also allows us to 
place tight constraints on the weakest sources; we find 
that <1% of sources are X-ray weak by a factor of >56 
(Aaox < -0.67). 

We have given examples of constraints on fractions of 
sources that are X-ray weaker than some specified value 
of Aaox ■ Now we wish to constrain this fraction for ar- 
bitrary levels of X-ray weakness. For the combined sam- 
ple of 217 SDSS and BQS non-BAL QSOs, wc calculated 
an upper limit on the fraction, fweaki^), of optically- 
selected QSOs that have Aaox < The distribution of 
fweak{x) is shown in Figure O for radio-quiet, non-BAL 
QSOs in B combined with the BQS QSO sample, and 
also for the B sample alone. 

To generate Figure [5l we used the following method. 
For each measured value of Aaox in our sample, we 
determined the number, N , of objects that had a smaller 



X-ray Emission from SDSS Non-BAL QSOs 



7 



value of Aaox- We used the binomial distribution to 
find an upper limit on the fraction, fweak, of such objects 
such that we would be 90% likely to observe at least A^+l 
objects that are equally X-ray weak or X-ray weaker if 
the true fraction of such objects is fweak- In Figure[5l we 
plot the upper limit on the fraction, fweakix), of objects 
that have Aaox < x. 

Figure [6] shows the same analysis performed for the 
sample of radio-quiet, non-BAL QSOs from B split into 
three subsamples according to their UV luminosities, 
LjgQQ^. We draw the tightest constraints from the lowest 
range of luminosities, although the constraints are highly 
influenced by differences in the relatively small number 
of sources with Aaox S —0.1. We do not find any clear 
evidence that the fraction of X-ray weak objects depends 
on luminosity. 

3.3. UV Absorption in Non-BAL QSOs 

BAL QSOs generally have reddened continua com- 
pared to those of non-BAL QSOs (e.g., iReichard et al.l 
120031 ). Apart from the presence of BALs, continuum 
reddening alone is no t known to correlate strongly with 
X-ray weakness, and iBrandt et al.l (|2000D found no sig- 
nificant difference between the continuum slopes of X-ray 
weak and non- X-r ay weak QSOs. U sing ROSAT X-Tay 
observations, Strat eva et al.l ()2005f ) found no trend of 
redshift-corrected g -~i colors with aox ■ We find no evi- 
dence in our sample for a correlation between Aaox and 
the ratio between UV flux densities at 1400 and 2500 A, 
which we use as a proxy for intrinsic reddening. Apart 
from reddening, this result also suggests that the ratio 
of UV/X-ray emission is not strongly dependent on the 
slope of the UV continuum (between 2500 and 1400 A) 
once the trend between UV and X-ray luminosities is 
factored out with Equation [H 

As described in §2.2) we have implemented an algo- 
rithm to search for narrow Mg II absorption doublets 
that would be associated with intervening systems or 
with absorbing material intrinsic to the QSO (or its host 
galaxy) . Figure [7] shows the distribution of Aaox for 
radio-quiet, non-BAL QSOs that have no Mg II dou- 
blets (top panel) compared to those that have at least 
one narrow Mg II doublet (bottom panel). We have only 
considered values for sources at redshift z < 2.28 (so 
that intrinsic Mg II at rest would be visible in the spec- 
trum) and in which the optical/UV continuum is > 4cr 
above zero for at least 50% of the spectral bins in or- 
der to minimize any bias due to the fact that sources 
with low signal-to-noise may have undetectable absorp- 
tion lines. The Gehan and logrank tests implemented in 
ASURV find no strong evidence that the distributions of 
Aaox differ between systems with and without narrow 
Mg II absorption lines. If X-ray absorption is related to 
the presence of narrow Mg II absorption doublets, the 
association is more subtle than we can detect with our 
current sample size, or may be driven by absorption fea- 
tures which are weaker than our EW threshold ( §2.2p . 

3.4. X-Ray Weakness and UV Emission Lines 

In order to investigate any relations between UV emis- 
sion line properties and X-ray weakness, we test for cor- 
relations between Aaox and emission line EW, FWHM, 
and central wavelength using Kendall's Tau test as imple- 
mented in ASURV. We consider only the C IV A1549 and 



C III] A1909 emission lines because other lines were not 
always present in the optical/UV spectra, and also were 
more difficult to constrain due to nearby spectral struc- 
ture. We define the emission line wavelength to be the 
center of our best-fit Voigt profile. We omit one source, 
J152156. 48 -1-520238. 4, from our fits because it has an un- 
usual emission line structure and is an outlier from our 
sample. The results of our correlation tests are shown 
in Table HI We report the results of correlation tests for 
radio-quiet sources including and excluding BAL QSOs 
in the sample. 

We see no strong (>99% confidence) correlations be- 
tween Aaox and the C IV FWHM, the C III] EW, or the 
C III] wavelength for radio-quiet sources. The C IV EW 
and Aaox are correlated at >99.99% confidence (Fig- 
ure El). The C III] FWHM and Aaox are also signifi- 
cantly anti-correlated at > 99.9% confidence (Figure [5]). 
The C IV wavelength is correlated with Aaox at 97-99% 
confidence (FigurefTO]). although this correlation weakens 
if we take the C IV wavelength with respect to the C III] 
wavelength rather than 1549.5 A in the rest frame deter- 
mined by the SDSS redshift. Overall, we find that the 
X-ray weakest sources (including BAL QSOs) have the 
weakest C IV emission, the broadest C III] emission lines, 
and blueshifted C IV emission. 

Fitting the trends that were found to have correlations 
at >99% confidence, we obtain: 

Aaox = (0.21 ± 0.03) log(C/V EW) + (-0.36 ± 0.0<j^) 



Aaox = (-0.15 ± 0.04) log(C///] EWHM) + (0.19 ± 0.06), a(i74 



Aaox = (0.004 ± 0.002)(C/V A - 1549.5) + (-0.026 ± 0.00®) 

where the EW, FWHM, and central wavelength are in 
Angstroms. 

SDSS redshift determination depends on QSO emission 
lines, and therefore the redshifts we use may have some 
dependence on the C IV emission wavelength. As both 
line wavelength and Aaox depend on redshift, correla- 
tions between the two can be affected by redshift deter- 
mination. The degree of change in Aaox due to changes 
in the redshift should be small compared to the trend in 
Figure [51 but we cannot fully account for the effects of 
emission line shifts on the SDSS algorithm. Aaox is also 
marginally correlated (at «96% confidence) with the dif- 
ference in the central wavelengths fit to C IV and C III] 
emission, and this correlation should be less sensitive 
to the estimated redshift. However, because of weaker 
confidences and redshift determination issues, we do not 
draw physical conclusions from the observed correlation 
between emission wavelength and Aaox ■ 

In order to illustrate more clearly the effect that corre- 
lations with line properties can have on Equation 31 we 
have repeated the fit described in [J33]with the C IV EW 
as a second independent variable. Our best fit is: 

aox{Lu, CIVEW) = (-0.158 ± 0.035) log(L25ooi) + (0-210 ± 0.04i; 

Equation [51 characterizes typical X-ray luminosities of 
optically-selected, radio-quiet, non-BAL QSOs (in the 
luminosity range of our SDSS sample) based on the UV 
luminosity and (relatively unknown) physical effects that 
determine the C IV emission EW. While the coefficients 
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on -Zj250oa ^ both non-zero at «5ct sig- the lower UV lumin osity (L 

nificance, the magnitude of the luminosity term is about 
fifteen times larger than the EW term for the luminosi- 
ties and line EWs in our sample. We also caution that 
the UV luminosity range in our sample is limited, and 
we cannot be sure that Equation [9] extends to other UV 
luminosities. 

3.5. Constraints on UV/X-ray Variation 

Even when radio-loud and BAL QSOs are removed 
from samples, there is still a significant amount of scatter 
about the best-fit UV/X-ray luminosity relation (Equa- 
tion[3]) . The scatter may be caused by additional physical 
factors, some of which have been considered in ii3.3l and 
§3.41 In this section, we place a limit on the maximal 
extent to which the UV and X-ray luminosities may vary 
with respect to each other. To determine this limit, we 
make the "worst-case" assumption that the scatter is en- 
tirely caused by source variability. We also assume for 
statistical purposes that Aaox is normally distributed, 
even though an Anderson-Darling test suggests this may 
not be perfectly accurate fi )3.ip . Figure S] shows that the 
true distribution does not differ too strongly from a nor- 
mal distribution, and this assumption should not greatly 
affect our result. Finally, our limit is also weakened 
by the non-simultaneity of our UV and X-ray measure- 
ments. Even if the UV and X-ray luminosities vary syn- 
chronously, non-simultaneous measurements will show 
increased scatter over simultaneous observations. For all 
these reasons, our analysis only gives an upper limit on 
the relative variation between UV and X-ray luminosi- 
ties. 

Following iMaccacaro et al.l (|1988f ) , we have minimized 
the likelihood function 

Aw4t ^ exp(-(Aaox. - t^fmaf + of^O) 

in order to determine the intrinsic mean (fi) and rms (cr) 
for our Aaox values, each having an estimated measure- 
ment error ctj. For radio-quiet, non-BAL QSOs in sample 
B, we find ^J. = 0.005 ± 0.015 and a = 0.115 ± O.OII. In 
fact, the square root of the unbiased sample variance is 
0.125, demonstrating that the influence of the measure- 
ment error is small compared to the spread of Aaox- 
The careful (unbiased) construction of sample B has al- 
lowed us to improve constraints on the spread of Aaox 
our cr is signif icantly sm aller than the aox rms of saO.2 
found by Steffen et al.l 12006.) for QSOs in their sample 
with 30 < log(L25ooA) < 

Assuming the width of the distribution of Aaox is 
purely due to luminosity variation (and measurement er- 
ror), we can use the intrinsic rms to constrain QSO X-ray 
fluxes to vary upward or downward with respect to UV 
fluxes by a factor of only about ioO ii5/0-3838 « 2 for 
sources within la of the mean, /i. For sources within 2a 
and 3a of the mean, this corresponds to factors of «4 and 
Ri8, respectively. This estimate only accounts for X-ray 
variation with respect to the UV continuum. UV fluxes 
may also be varying along with the X-ray fluxes, so that 
the absolute X-ray variation is larger. Compared to this 
upper limit. X-ra y variation by a facto r of 30, such as that 
seen m Mrk 335 (|Grupe et ai.ll2007a^ ■ is rather unusual, 
unless the UV continuum has also varied strongly. How- 
ever, we note that our SDSS sample does not extend to 



, , . - 2500 A - 10'' erg s-i Hz-i, 

iGrupe et al.l l2007a[ ) of this source, and that variability 
properties may differ at lower luminosities. 

We have also performed this analysis calculating 
Aaox using Equation [9] in order to remove scatter in 
the UV/X-ray luminosity relation that is correlated with 
the C IV emission EW. These effects are small, how- 
ever, compared to the overall span of UV and X-ray lu- 
minosities, and the end result is only a slightly tighter 
constraint on variation. We obtain in this case /i = 
0.009 ± 0.014 and a = 0.103 ± 0.011. 

Nearby Seyfert 1 AGN typically show fractional vari- 
ability in their X-ray l uminosities by up to » 40% in 
X-rays over a few years ()Markowitz et al.]|2003( ). Apart 
from this general study, dramatic variation has been ob- 
served in individu al sources such as NGC 4051 (e.g., 
iLamer et al.l |2003[ ) . The variability characteristics of 
higher-luminosity AGN are not well-understood; recent 
studies have found some evidence that QSOs vary less at 
higher lumino sities, but tha t variability increase s at z « 
2 (e.g., .Manners et all 120021 : iPaolillo et aIll2004D. In in- 
dividu al cases, variation can be extreme. Ih emmer et al.l 
(|2005f l observed luminosity variation by factors of «2-4 
over 3-5 (rest frame) months in several QSOs at z > 4. 

In order to test whether variability could account for 
the spread of Aa px, we assume a UV variability am- 
plitudc of «30% (iGiveon et all [19991 : iKaspi et al.ll200l 
that is uncorrelated to the X-ray variation (~40%). We 
caution that these are simply estimates and that the vari- 
ation properties are poorly constrained. In some individ- 
ual cases, high-luminosity QSOs have shown variation 
of up to 70% in t he UV over 2-3 yr in the rest-frame 
(jKaspi et al.ll2007t ). and in X-rays by a factor of fti4 over 
73 days in the rest-frame (jShemmer et al.l[2005f ). Using 
these estimates, the width of the distribution of Aaox 
is largely (70-100%) attributable to intrinsic variation. 

4. DISCUSSION 

In this section, we discuss additional implications of 
the results presented in §??, including: the rarity of ex- 
tremely X-ray faint, optically-selected QSOs (t J4.ip : the 
relation of the anomalously X-ray faint QSO PHL 1811 
to our SDSS sample ( tj4.2p : and the dependence of UV 
emission line strengths on X-ray luminosity f t|4.3p . 

4.1. X-Ray Weakness in Optically- Selected QSOs 

Using X-ray data from the Einstein Observatory, 
lAvni fc TananbaumI (|1986[ ) derived an upper limit of 8% 
(at 95% confidence) for the fraction of optically-selected 
QSOs that could be so X-ray faint as to be a separate 
population from "ordinary," X-ray-emitting QSOs. At 
the time of their study, it was not known that BAL 
QSOs were anomalously X-ray weak. However, Figure [3] 
demonstrates that BAL QSOs generally live on the tail 
of the overall distribution of relative X-ray brightness, 
and would likely n ot be considered a disjoint population 
by the criteria of lAvni fc TananbaumI (|1986D . 

The high sensitivity of Chandra has enabled us to con- 
strain the fraction of QSOs that lie in the X-ray-weak 
tail of the distribution of "ordinary" QSOs. If a separate 
population of extremely X-ray weak QSOs exists, as in 
the hypothesis tested bv lAvni fc TananbaumI (|1986D . our 
results also provide upper limits on the fraction of such 
sources in optically-selected samples. Although we can- 
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not rule out the possibility that a small fraction of QSOs 
are extremely intrinsically X-ray faint, we have deter- 
mined that such sources are at least very unusual in the 
SDSS. 

In fact, we find evidence that few (<2%) SDSS QSOs 
are relatively X-ray weak by even a factor of 10. When 
BQS QSOs are included in the sample, we find that <1% 
of sources are X-ray weak by a factor of <56 {Aaox < 
—0.67). Thus, if accretion disk coronae are responsible 
for X-ray emission in QSOs, they appear to be present 
and effective in most, if not all, luminous AGN. Figure[5] 
shows our estimates of upper limits on the fraction of 
radio-quiet, non-BAL QSOs that are X-ray weaker than 
a given value of Aaox- These limits are driven by the 
sizes of our samples. 

Figure [5] allows the calibration of wide, shallow X-ray 
surveys to constrain the fraction of optically-selected 
QSOs that should be detected. The current upper lim- 
its indicate that at most a few percent of optically- 
selected radio-quiet, non-BAL QSOs will be missed at 
Aaox < —0.4. We have achieved a 100% detection frac- 
tion for the sources in sample B with exposure times 
as short as 2.5 ks. Future shallow, wide X-ray surveys 
should therefore be able to detect nearly all non-BAL 
QSOs with UV luminosities as faint as those in the SDSS 
in relatively short exposure times. 

4.2. How Unusual Is PHL 1811? 

A multiwavelength study of the nearby [z = 0.192) 
narrow- line QSO PHL 1811 has revealed that this 
source has both unus ual X-ray and optical/UV spectra 
(|Leighlv et al.l l2007bl la[). It has Aaox ~ -0.7, corre- 
sponding to X-ray weakness by a factor of almost 70 
compared to QSOs with the same UV luminosity. There 
is little, if any, intrinsic X-ray absorption, the UV spec- 
trum shows no broad absorption lines, and broad opti- 
cal/UV emission lines are very weak. C IV A1549 has 
an EW of only 6.6 A, and only a limit of <1.7 A can be 
placed on the C III] EW. 

If PHL 1811 were at a redshift z = 2, its contin- 
uum properties would justify its inclusion in our sam- 
ple. Its blue color would identify it as a QSO candidate 
for SDSS spectrosc o py, acc ording to the criteria s pecified 
by [Richards et"aLl (|2002bO . and iLeighlv et all i|2007aD 
note (in their §4.5) that it would have been identified 
as a PG QSO, if it had been in that survey's footprint. 
The 2500 A monochromatic lumin osity of PHL 1811 is 
-^250oi ~ 10^° erg s-i Hz"! (jLeighlv et al.ll2007bl 
§4), which is intermediate for our sources (Figure [5]). 
Thus, our survey can be used to constrain the frequency 
of intrinsically X-ray weak objects like PHL 1811 in the 
Universe. From Figure O it appears that the weak in- 
trinsic X-ray emission of PHL 1811 makes it an unusual 
source that is not representative of a significant frac- 
tion of optically-selected QSOs. If the X-ray weakness of 
PHL 1811 is due to an abnormal a c cretion disk corona, as 
discussed in §5.1 of ILeighlv et al.l (|2007bl ). such physical 
states are apparently rare in optically-selected AGN. 

The weak, blueshifted C IV A1549 emission observed in 
PHL 1811 is qualitatively consistent with the trends we 
observe for X-ray weaker objects f §3.4p . However, it is 
still an outlier compared to the emission line properties 
of our sample. For example. Equation [11] predicts a C IV 
emission EW of 17 A based on UV and X-ray luminosi- 



ties, while ILeighlv et al] (|2007af ) actually measure 6.6 A. 
If the X-ray luminosity were to increase, this would also 
increase the predicted emission line strength, making the 
discrepancy worse. PHL 1811 is therefore an outlier from 
our sample according to both continuum and emission- 
line properties, for physical reasons that are currently 
unknown. 

4.3. Broad Emission Lines Correlated with X-Ray 
Brightness 

While the strong UV absorption in BAL QSOs is 
known to be associated with X-ray weakness, it is not 
clear how strongly other properties evident in the op- 
tical/UV spectrum are associated with X-ray weakness. 
We do not find any strong correlations between X-ray 
weakness and UV reddening or the presence of narrow 
absorption line systems. We do find significant correla- 
tions between Aaox and the C IV emission EW, the 
C III] emission FWHM, and the wavelength shift of the 
C IV emission line. 

It has been known for some time that C IV emis- 
sion EWs are anti-correlated with the UV c ontinuum 
strength. This is the well-known Baldwin effect (jBaldwinI 

Il977f) . typically parameterized as EW oc Lfjy, where 
Ljjv is the UV luminosity. For an ensemble of QSOs, the 
"global Bal dwin effect" has a logarithmic slope f3 « —0.2 
(IKinnev et al. 1990). For a sample of 105 SDSS QSOs, 
IWilh ite et al.i (i200Ql find P = -0.22 ± 0.03. For our 
sample, we find /3 = —0.24 ±0.03, a bit steeper than but 
consistent with previous s tudies. The co mparison is com- 
plicated by the fact that IKinnev et al.l (fl99Ql measured 
the UV luminosity at 1550 A rather than at 2500 A, but 
comparison of these luminosities for our data indicates 
that the parameter f] should not be greatly affected. 

Equation |S| indicates that the C IV emission EW may 
also be related to X-ray luminosities. This is not unex- 
pected in a photionization-driven scenario, as the K and 
L shells of C IV have ionization potentials of «320 and 
«64 eV, respectively. Fitting the emission EW as a func- 
tion of both UV and X-ray luminosity for sample B, we 
find: 

\og{CIVEW) = (-0.41 ± 0.07) log(i25ooi) + (0-29 ± 0.06) log(L2 ke 
According to Equation 1111 when X-ray luminosities are 
considered, the dependence on UV luminosities steepens 
somewhat to (3uv ~ —0.41. The coefficient on the X-ray 
luminosity is also significant, with Px ~ 0.29. The two 
coefficients have opposite signs, so that C IV emission 
weakens with increasing UV luminosity, but increases 
with X-ray luminosity. 

By definition, aox increases as L2 koV increases or 
-^250ol decreases. The opposite signs and roughly com- 
parable magnitudes of the coefficients f3x and Puv in 
Equation [Tl] indicate that the C IV emi s sion s trength 
is correlated with aox- In fact, iGreenI (Il998h previ- 
ously observed that the C IV emission EW depended 
more strongly on aox than on UV luminosity alone, 
and suggested that the emission line strength may there- 
fore be driven by the continuum shape (e.g., the slope 
between UV and X-ray luminosities characterized by 
ctox)- A correlation between C IV EW and the ra- 
tio of monochromatic lu minosities at 1 keV and 1350 A 
has also been found by IWang et al.l (|1998) . Such ob- 
servations correspond well with photoionziation models 
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of equilibrated gas, in which the shape of the ionizing 
continuum is an important parameter for determining 
ionization fractions. In fact, it has been argued that the 
Bald win effect may be driv en by the photoionizing flux 
(e.g.. iKorista fc Go ad 2004*). although further modeling 
is needed to determine whether the X-ray dependence 
of the Baldwin effect is consistent with photoionization 
scenarios. 

On the other hand, C IV emission line properties are 
know n to be related to t he Eddington ratio, L/LecIcI 
(e.g.. iBaskin fc Laoil [20051) . The Hp line is redshifted 
out of the SDSS spectral range, so we are unable to de- 
termine Eddington ratios for our sources using H/s line 
diagnostics. As an alternative to photoionzation scenar- 
ios, the relation between C IV EW and UV/X-ray lumi- 
nosities may be due to an underlying physical parame- 
ter, such as Eddington ratio or perhaps orie ntation (e.g., 
[Richards et"aIll2002aHGallagher et al.ll2005D . that drives 
the C IV emission and both the UV and X-ray luminosi- 
ties. It has also been suggested that C IV blueshifts, 
the Baldwin effect, and relative X-ray brightness have a 
common physical origin, and that the underlying phys- 
ical parameter governs a spectrum of QSO properties. 
BAL QSOs wo uld fall on the weak-emission extreme of 
this spectrum (lRichardsll2006l and references therein). 

5. CONCLUSIONS AND FUTURE WORK 

We have analyzed the relationship between X-ray lumi- 
nosity and UV properties of SDSS DR5 QSOs at redshift 
1.7 < z < 2.7. At these redshifts, we are able to identify 
C IV BAL QSOs in our sample. This has enabled us to 
place strong constraints on the fraction of intrinsically 
X-ray weak QSOs, as well as to search for additional 
trends beyond the well-studied relationship between UV 
and X-ray luminosities. In particular, we find that: 

1. Non-BAL, radio-quiet QSOs which are very (in- 
trinsically) X-ray weak are rare in optically- 
selected samples, with <2% of SDSS QSOs hav- 
ing Aaox < -0.4 and <1% of SDSS/BQS QSOs 
having Aaox < —0.67. Figure [5] provides upper 
limits on the fraction of sources that may be rela- 
tively X-ray weaker than a given value of Aaox ■ 

2. The rms of the Aaox distribution for radio-quiet, 
non-BAL QSOs is about 0.1, corresponding to a 
factor of f»2 spread in X-rays relative to that esti- 
mated from the UV luminosity. This places an up- 
per limit on typical X-ray variability with respect 
to the UV continuum. 

3. While the amplitude of and relation between 
UV and X-ray variation in QSOs is not well- 
understood, estimates based on recent studies in- 
dicate that most, if not all, of the observed spread 
of Aaox can be attributed to variability. 
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4. The distribution of Aaox may not be normally 
distributed (at 97.5% confidence) in our sample B. 
Perhaps physical processes such as X-ray absorp- 
tion are significant in some sources, broadening the 
wings of the distribution. 

5. We find no strong evidence that reddening or the 
presence of narrow Mg II absorption is related to 
X-ray weakness for non-BAL QSOs. 

6. UV emission line properties such as C IV EW, 
C III] FWHM, and perhaps C IV wavelength are 
correlated to relative X-ray brightness. 

7. The C IV emission EW depends on both UV and 
X-ray luminosity. The physics that drives the 
global Baldwin effect is apparently associated with 
X-ray emission as well as UV emission. 

8. Even after correcting for secondary trends (such as 
weak C IV emission, which is associated with rela- 
tive X-ray weakness), objects that are as intrinsi- 
cally X-ray faint as PHL 1811 are rare. 

We have quantitatively shown that luminous X-ray 
emission is essentially a universal property of optically- 
selected QSOs. Future studies can test whether this re- 
sult holds true for QSOs selected in other wavebands, 
such as the radio and infrared. New optical surveys will 
extend the luminosity range of our sample, while UV 
(and infrared) spectroscopy will allow BAL QSOs to be 
identified (and ^2500 A measured) for sources in 

a wider redshift range. The high detection rate for rel- 
atively short exposures in our sample has also demon- 
strated that wide, shallow X-ray surveys at high angular 
resolution are an effective way to study the X-ray prop- 
erties of bright, optically-selected QSOs. 
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TABLE 1 

Source Information^ 



SDSS 


z 


Obsld*^ 


X-Ray 


X-Ray Counts 


X-Ray 




Sample 




log(i?*) 


CIV BIq 


Narrow Mgll 


log(L, 


Iog{L,-,,„ 


log(L2 kcv)'^ 




Source 






Exposure (ks)"^ 


Soft/Hard'* 


Detected?'^ 










(km s-l) 


Systems^ 






J000654.10 - 001533.4 


1.73 


4096 


4.5 


34/10 




1 







< 0.19 





3 


31.01 


31.23 


26.89 


-1.67(0.03) 


J000659.28 - 001740.8 


2.02 


4096 


4.5 


19/5 




1 




1 


< 0.75 








30.70 


30.80 


26.77 


-1.55(0.04) 


J000717.73 - 002811.3 


1.89 


4096 


4.5 


11/0 




1 







2.18 








30.03 


30.31 


27.12 


-1.22(0.07) 


J001130.55 + 005550.7 


2.31 


0403760301 


25.4 


34/29 




1 







< 0.28 


928 





31.02 


31.37 


26.33 


-1.94(0.03) 


J001247.12 + 001239.5 


2.15 


4829 


6.7 


9/3 




1 




1 


< 0.54 








30.82 


31.05 


26.61 


-1.71(0.06) 


J001306.15 + 000431.9 


2.16 


4829 


6.7 


13/3 




1 







< 0.37 
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31.17 


31.23 


26.44 


-1.84(0.06) 


J002025.22 + 154054.6 


2.01 
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19.9 
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3.20 
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-1.60(0.01) 
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1.76 
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^ The full version of this table is available in the electronic edition online.^ The Chandra or X M M- Newton observation identification number. Ten-digit numbers correspond to XMM-Newton observations, while short 
correspond to Chandra observations.'^ The effective X-ray exposure reported by the CIAO or SAS toolchains for the extraction of the source region.*^ The observed-frame soft (0.5—2 keV) and hard (2-8 keV) counts in 



region.^ A "1" indicates an X-ray detection; "0" indicates no detection.^ A "1" indicates the source is in sample B; "0" indicates it is not in sample B.^ The number of narrow Mg II absorbtion systems found using i 
described in ^2.2| ^ Monochromatic luminosities are given in units of erg s~^ Hz~^.' The error, assumed to be dominated by the error on X-ray luminosity, is reported in parentheses. 
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UV Emission Line Measurements'' 
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J002155.52 + 001434.3 


19.5 


1548.7 






15.9 


1905.0 


34.0 


J002308.74 + 002239.7 


34.5 


1545.4 




16.1 


10.0 


1902.4 


30.5 


J002331.21 - 011045.6 


20.8 


1541.0 




17 9 


18.2 


1896.1 


42.1 


.1002825.59 + 003500.1 


24.7 


1552.1 




25.0 


30.4 


1915.4 


49.4 


J002917.31 - 002540.9 


80.1 


1547.7 




17.7 


36.2 


1907.9 


22.4 


J002954.96 - 001053.4 


13.7 


1544.8 




28.8 


12.9 


1901.0 


30.6 


J003131.44 + 003420.2 


59.6 


1545.6 




32.5 


19.4 


1901.2 


47.7 


J003135.56 + 003421.3 


46.9 


1537.1 




22.1 


36.3 


1895.4 


34.0 


J003922.44 + 005951.7 


31.3 


1551.4 




13.0 


6.3 


1912.5 


23.9 


J004206.18 - 091255.7 


32.5 


1546.6 




22.3 


22.2 


1904.8 


29.5 


J004349.50 + 003930.2 


87.0 


1545.6 




33.2 


41.8 


1898.3 


55.8 


J004526.26 + 143643.5 


39.9 


1544.1 




29.0 


17.0 


1906.0 


41.8 


J004527.68 + 143816.1 


14.3 


1547.7 




19.9 


19.0 


1896.4 


50.5 


J005018.84 - 005438.0 


14.4 


1545.2 




31.2 


14.5 


1906.4 


41.0 


J005025.13 - 005718.3 


42.4 


1555.2 




21.7 


28.4 


1918.7 


32.6 


J005102.42 - 010244.3 


29.5 


1546.3 




36.2 


21.5 


1909.1 


43.5 


J005355.15 - 000309.3 


23.7 


1539.2 




27.6 


18.4 


1899.8 


39.0 



^ The full version of this table is available in the electronic edition online. 
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TABLE 3 
Sample Information 



Sample 


Total 


Radio-Quiet 


Radio- Loud 


Radio-Quiet 


Radio-Loud 




Size 


Non-BAL QSOs'' 


Non-BAL QSOs"^ 


BAL QSOs'' 


BAL QSOs'^ 


A 


536 


433/354 


34/32 


64/43 


5/5 


B 


163 


139/139 


10/10 


13/11 


1/1 



^ The number of sources is reported as "a/6," where a is the total number of sources of the 
specified type, and 6 is the number of those sources which were detected in X-rays. 



TABLE 4 

Tests for Emission Line Correlations with Aaox 





Relation 


Confidence 
BAL QSOs Included 


Confidence 
Only Non-BAL QSOs 


C IV EW 


>99.99 


>99.99 


C IV Wavelength 


99.2 


97.0 


C IV FWHM 


67.5 


32.0 


cm 


EW 


94.8 


85.2 


cm 


Wavelength 


31.8 


36.2 


cm 


FWHM 


>99.99 


99.96 


CIV 


- C III] Wavelength 


90.3 


95.8 




z 

Fig. 1. — The redsliifts and absolute i magnitudes, Mj, for our 536 sample A sources at 1.7 < z < 2.7 (filled circles) compared to the 
entire SDSS DR5 QSO catalog (gray points). We also indicate (with filled squares) the BQS QSOs at z < 0.5 which were included in the 
SDSS DR5 QSO catalog. 



X-ray Emission from SDSS Non-BAL QSOs 



15 




log 



10^ 



n 



erg 
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Fig. 2. — O-ox measured for all radio-quiet, non-BAL QSOs plotted against the monochromatic luminosity at 2500 A. Radio-quiet, non- 
BAL QSOs in the subsample B are all detected and are plotted with open squares, while those detected in A (but not B) are plotted with 
filled squares. In cases where sources in A were not detected in X-rays, we have plotted upper limits, indicated with arrows. The dashed 
line shows the fit that Just et al. (2007) obtain ed fr om their sample of objects, while the s olid line shows the bes t fit to the radio-quiet, 
non-BAL QSOs in our sample B, described in i|3.1l A cross marks the values measured bv lLeighlv et all l)2007bl ) for the unusual source 
PHL 1811 (described in §??). PHL 1811 is not included in our SDSS sample, but would be an outlier if it were. 




O 



-0.5 







RL, BA[ 



0.5 



^a 



ox 



Fig. 3. — Aoox for the sources in our sample. From top to bottom, the panels show the distribution of Aaox for radio-quiet sources 
with no BALs, radio-loud sources with no BALs, radio-quiet sources with BALs, and radio-loud sources with BALs. 
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RQ, No BAL 
Sample B 




-0.5 



Fig. 4. — Aap x for the radio-quiet sources with no BALs in sample B. The solid Une shows the best Gaussian fit to the distribution, 
discussed in i|3.1l The x-axis is chosen to match that of Figure [5] 




X 

Fig. 5. — Upper limits on the fraction fu,eak{^) of sources which have Aaox < 2;, calculated as described in i|3.2l The dashed line 
connecting circles shows results for sample B combined with BQS QSOs. The dotted line connecting crosses shows results for sample B 
alone. The j;-axis is logarithmic. 




X 

Fig. 6. — Similar to Figure [5] we plot upper limits on the fraction f^eaki^) of sources which have Aapx ^ We use radio-quiet, 

non-BAL QSOs from sample B, and split the combined sample into three subsamples according to i2500A' '^^^ shows the plot symbols 
that correspond to each UV luminosity range. The y-axis is logarithmic. 
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Fig. 7. — The distribution of Aaox for radio-quiet QSOs with no BAL and no identified narrow-line Mg II absorbers (top panel) and 
for sources with at least one narrow-line absorber {bottom panel). 
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Fig. 8. — Aapx plotted against rest-frame EWs for C IV emission lines. All detected, non-BAL, radio-quiet QSOs are plotted as filled 
squares, with upper limits given by arrows. Detected, radio-quiet BAL QSOs are plotted as open circles, with upper limits for non-detected 
radio-quiet BAL QSOs designated by arrows inside the circles. The plot ranges clip off a small number of points. The solid line indicates 
the best fit to the data. 
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Fig. 9.— Same as FigureE] but showing the FWHM of C III] A1908 on the x-axis. 
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Fig. 10. — Same as Figure|S] but showing the central wavelength of C IV A1549 on the x-axis. 



